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ABSTRACT 


Using  archival  data  from  the  MTT  X-Ray  Observatory  on  board  the  SAS- 
3  satellite,  we  have  produced  refined  X-ray  light  curves  of  the  X-ray  nova 
A0620-00  during  its  initial  outburst  and  subsequent  decay.  We  use  data  from 
similar  observations  of  the  Crab  Nebula  to  generate  elevation  dependent 
response  functions  for  the  SAS-3  detectors  and  to  develop  objective  criteria 
for  perfecting  the  A0620  data  in  the  energy  ranges  1.3-5  keV  and  5-13 
keV.  We  then  analyzed  the  A0620  light  curves  and  the  hardness  ratio  for 
evidence  of  periodicity.  In  particular,  we  searched  for  evidence  of  the  7.75h 
orbital  period  found  in  optical  observations  by  McClintock  and  Remillard 
(1986).  Using  a  Monte  Carlo  simulation,  we  showed  that  the  3a  upper  limit 
on  the  amplitude  of  oscillation  in  the  hardness  ratio  is  3.2%.  We  also 
searched  for  the  7.8d  period  reported  by  Matilsky  et  al.  (1976).  We  found 
three  cycles  of  a  7.8d  oscillation  in  a  22  day  section  of  the  data  with  a  ~  40% 
amplitude.  However  in  the  remainder  of  the  data,  in  excess  of  200  days,  we 
were  able  to  place  3a  upper  limits  on  the  intensity  amplitude  of  2.2%  before 
the  section  and  7.3%  after.  We  therefore  find  no  evidence  to  substantiate  a 
persistent  7.8d  period.  We  conclude  that  the  system  exhibits  little  to  no 
persistent  X-ray  periodicity  at  any  frequency.  We  also  conclude  that  high 
inclination  angle  models  of  the  system  where  the  accretion  disk  is  fed  via 
Roche  lobe  overflow  are  not  appropriate. 
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Chapter  1:  Introduction 


On  August  3, 1975  the  British  satellite  Ariel  V  discovered  the  X-ray  nova 
A0620  (Elvis  et  al.  1975).  The  MIT  X-Ray  Observatory  on  board  NASA’s 
SAS-3  satellite  began  observing  it  on  August  8,  1975  and  continued  to 
observe  it  fairly  regularly  until  February  13, 1976.  The  intensity  of  the  nova 
rapidly  rose  to  a  peak  value  and  for  two  months  it  was  the  brightest  object  in 
the  X-ray  sky.  Using  the  X-ray  position  determined  by  Doxsey  et  al.  (1976), 
Boley  and  Wolfson  (1975)  identified  the  optical  counterpart  and  reported 
that  the  X-ray  nova  was  accompanied  by  an  outburst  in  the  visible  spectrum: 
the  optical  counterpart  brightened  from  18th  to  12th  magnitude.  Following 
the  initial  discovery  there  were  a  flurry  of  observations  across  the 
electromagnetic  spectrum.  Craft  (1975),  Scott  (1975),  and  Owen  et  al. 
(1976),  identified  a  radio  counterpart.  The  nova  was  observed  in  the  radio 
by  Little  et  al.  (1976).  It  was  observed  in  the  infrared  by  Citterio  et  al. 
(1976)  and  Kleinmann  et  al.  (1976).  Carpenter  et  al.  (1976),  Lloyd  et  al. 
(1977),  and  Tsumeni  et  al.  (1977)  observed  it  in  the  optical.  The  Ariel  V  X- 
ray  spectrum  was  observed  by  Ricketts  (1975)  and  Kaluzienski  et  al.  (1977) 
discussed  the  decay  of  the  X-ray  intensity  as  observed  by  Ariel  V.  The  SAS-3 
X-ray  observations  were  analyzed  by  Matilsky  et  al.  (1976),  MB  hereafter. 
Many  spectroscopic  studies  were  done  including  those  by  Coe  et  at.  (1976), 
Wheelan  et  al.  (1977),  Oke  et  al.  (1977),  and  Long  et  al.  (1978). 

After  15  months  the  system,  a  low  mass  binary  with  a  compact 
component,  returned  to  its  pre-outburst  brightness  and  has  remained 
quiescent  ever  since.  Using  the  Harvard  plates,  Eachus  et  al.  (1976)  found  an 
outburst  of  the  optical  counterpart  in  1917,  and  it  seems  likely  that  it  too  was 
accompanied  by  an  X-ray  nova.  Photometric  and  spectroscopic  studies  by 
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McClintock  and  Remillard  (1986),  MR  hereafter,  demonstrated  the  system 
has  a  high  mass  function  and  placed  a  lower  limit  on  the  mass  of  the  compact 
component  which  exceeds  the  theoretical  maximum  mass  of  a  neutron  star. 
This  study  sparked  interest  in  A0620-00  as  a  black  hole  binary  system. 

In  this  paper,  we  analyze  the  X-ray  data  obtained  from  the  detectors  on 
the  SAS-3  X-Ray  Observatory,  the  instrumentation  for  which  was  designed 
and  built  at  the  MIT  Center  for  Space  Research.  We  use  these  data  to  search 
for  evidence  of  the  7.75h  orbital  period  in  the  X-ray  phenomena  and  to 
reexamine  the  evidence  for  the  7.8d  period  reported  by  MB.  In  Chapter  1, 
we  present  background  information  on  black  holes  and  accretion  disk 
systems  relevant  to  a  study  of  A0620  as  a  black  hole  candidate.  We  also 
review  the  three  major  black  hole  candidates  and  recount  various 
researcher's  observational  evidence  in  favor  of  the  candidates  being  black 
holes.  In  Chapter  2,  we  discuss  how  we  processed  and  refined  the  SAS-3  data 
and  describe  how  we  used  data  from  Crab  Nebula  observations  to  assess  and 
correct  systematic  defects  in  the  processing.  In  Chapter  3,  we  apply  the 
criteria  and  procedures  of  Chapter  2  to  A0620  to  generate  refined  light 
curves  that  allow  accurate  detrending  and  thereby  a  more  sensitive  search  for 
periodicities.  In  Chapter  4,  we  discuss  the  theory  behind  the  period 
searching  techniques  employed  and  apply  these  techniques  to  generate  limits 
on  the  observed  X-ray  periodicities  of  A0620. 


1.1  Black  Holes 

When  a  star  exhausts  its  available  nuclear  fuel,  thermal  pressure  cannot 
stop  its  contraction.  However,  electron  degeneracy  pressure  can  sustain  an 
evolved  star  against  gravitational  collapse  as  long  as  its  mass  does  not  exceed 
the  Chandrasekhar  limit  of  1.4  Mq.  Stars  supported  solely  by  electron 
degeneracy  pressure  are  known  as  white  dwarfs.  They  typically  have  a  mass 
of  about  one  Mq  and  a  radius  of  about  5000  km. 

If  a  white  dwarfs  mass  exceeds  the  Chandrasekhar  limit,  the  electron 
degeneracy  pressure  is  insufficient  to  resist  the  gravitational  attraction.  The 
star  collapses.  Electrons  and  protons  merge  via  inverse  beta  decay  and 
become  neutrons.  Total  collapse  is  prevented  by  neutron  degeneracy 
pressure.  Stars  supported  by  neutron  degeneracy  pressure  are  called  neutron 
stars.  They  are  observed  as  radio  pulsars  and  as  the  X-ray  emitting 
components  of  binary  X-ray  sources. 

The  limits  of  stability  of  neutron  stars  are  not  completely  understood. 
Rotation,  and  the  equation  of  state  (which  is  not  precisely  known)  are 
essential  determinants  of  the  critical  mass.  However,  within  the  framework 
of  general  relativity  one  can  derive  an  upper  limit  based  on  considerations  of 
causality,  which  restricts  the  equation  of  state  so  as  to  yield  a  speed  of  sound 

less  than  the  speed  of  light.  According  to  Shapiro  and  Teukolsky  (1983)  (ST 
hereafter),  the  mass  limit  is  probably  about  2.7  Mq  with  an  absolute  causality 

limit  of  ~  3  -  5  Mq. 


4  Chapter  1:  Introduction 

1.1. a  Characteristics 

In  theory  four  observables  characterize  a  black  hole:  charge,  angular 
momentum,  radius,  and  mass. 

Since  celestial  objects  are  embedded  in  conducting  plasmas,  they  cannot 
sustain  a  condition  of  substantial  electrical  charge  under  any  circumstances. 

The  angular  momentum  of  a  black  hole  affects  the  position  of  the  event 
horizon  and  can  in  principle  create  certain  interesting  mechanical  and  optical 
effects.  However,  these  effects  are  difficult  if  not  impossible  to  observe  at 
astronomical  distances. 

The  radius  of  a  black  hole  is  its  event  horizon  at  which  unique  optical 
phenomena  could,  in  principle,  be  detected  by  visual  means.  However  the 
angular  resolution  required  for  the  nearest  possible  black  hole  is  hopelessly 
unattainable.  It  appears,  therefore,  that  the  only  feasible  method  for 
detecting  and  characterizing  a  black  hole  at  astronomical  distances  is  through 
measurement  of  its  mass  and  its  identification  as  a  collapsed  star. 

If  the  object  happens  to  be  in  a  binary  system,  then  we  can  observe  its 
effect  on  its  visible  companion.  Using  Kepler’s  third  law,  one  can  relate  the 
companion's  semimajor  axis,  period,  and  mass  (obtainable  through  its 
spectral  type)  to  the  mass  of  the  compact  object.  A  convenient  formula  for 
relating  the  observed  orbital  parameters  to  the  compact  object’s  mass  is  the 
mass  function 


f(M„  M*  i)  = 


(M2  sin  i)' 


(Mj  +  Mj) 


Pv? 

2ttG 


(l.l.a.l) 


1.1  Black  Holes 


5 


where  Mi  is  the  mass  of  the  optical  companion  star,  M2  is  the  mass  of  the 
compact  object,  i  is  the  inclination  angle,  and  vi  is  the  projection  of  the 
orbital  velocity  along  the  line  of  sight  of  the  observer  (see  figure  1.1).  The 
mass  function  depends  only  on  the  observables  P,  vi,  and  i  and  is  equal  to  the 
minimum  possible  value  of  M2  which  occurs  when  Mi  =  0  and  sin  i  is  a 
maximum  (i  =  90°).  In  chapter  1 .4  I  will  show  how  others  have  applied  the 
mass  function  mass  determination  method  to  estimate  the  masses  of  CYG  X- 
1,  LMC  X-3,  and  A0620. 


ai  % 

■< - >- 

a 


Figure  1.1:  The  parameters  of  the  binary  black  hole  system.  Mj  is  the  star,  M  is  the 
hole,  and  i  is  the  inclination  angle  with  respect  to  the  observer. 
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1.1. b  Useful  Results  from  General 
Relativity 

From  the  minimization  of  the  general  relativistic  action,  it  can  be  shown 
(ST)  that  the  specific  energy,  e,  of  a  test  particle  in  a  circular  orbit  about  a 
non-rotating  black  hole  is  given  by 


r  -  2M 
Vr(r  -  3M) 


(l.l.b.l) 


where  r  is  the  radius  of  the  orbit,  M  is  the  mass  of  the  central  body  and  the 
"geometrized"  units  G  =  c  =  1  have  been  used.  While  circular  orbits  exist  all 
the  way  down  to  r  =  3M,  they  are  only  stable  down  to  r  =  6M.  Within  r  =  6M 
a  small  perturbation  will  cause  the  particle  to  spiral  rapidly  into  the  central 
object.  Substituting  r  =  6M  into  equation  (l.l.b.l)  one  obtains  the  binding 
energy  of  the  last  stable  orbit 


Ebinding  =  (m  -  E)/m  =  1  -  e  =  5.72%  (1.1  .b.2) 

This  represents  the  fraction  of  its  rest  mass  energy  that  a  particle  originally  at 
rest  at  infinity  releases  as  it  spirals  down  to  the  last  stable  circular  orbit  and 
before  it  plunges  into  the  hole.  The  release  of  gravitational  energy  is 
evidently  a  very  efficient  means  of  energy  production,  a  conclusion  that 
motivates  the  idea  that  accretion  by  a  black  hole  may  explain  the  huge 
luminosities  that  some  compact  objects  generate. 

More  complicated  is  the  Kerr  (rotating)  black  hole.  We  define  the 
angular  momentum  per  unit  mass  a  =  J/M.  When  a  =  0  one  returns  to  the 


1.1  Black  Holes 
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Schwarzschild  (non-rotating)  case.  If  a  =  M  then  the  black  hole  is  said  to  be 
maximally  rotating.  As  in  the  Schwarzschild  case,  one  can  solve  for  the 
specific  energy  which  can  be  shown  to  be 


r^  -  2Mr  ±  aV  Mr 
rV r2  -3Mr  ±  2a V Mr 


(l.l.b.3) 


where  the  upper  sign  is  for  a  corotating  (direct)  orbit  and  the  lower  sign  is 
for  a  counterrotating  (retrograde)  orbit.  A  quantity  of  particular  interest  is 
the  binding  energy  of  the  last  stable  circular  orbit  for  a  maximally  rotating 
black  hole.  When  a  =  M,  the  radius  of  the  marginally  stable  orbit  is  rms  =  M 
(direct),  or  9M  (retrograde).  For  a  direct  orbit  at  this  radius  and  angular 
momentum,  e  =  (1/3 )1/2  and  the  binding  energy,  which  as  before  represents 
the  fraction  of  rest  energy  released  by  a  particle  initially  at  rest  at  infinity,  is 


^binding  —  1  -  e  —  42.3%! 


(l.l.b.4) 


One  can  see  that  the  release  of  gravitational  binding  energy  is  a  phenomenally 
efficient  means  of  energy  production.  This  is  why  black  holes  have  been 
invoked  to  explain  objects  from  Quasars  to  active  galactic  nuclei. 

In  section  1.2  I  will  present  the  results  of  the  Keplerian  approximation  to 
an  accretion  disk  about  a  black  hole  and  show  that  the  disk  strongly  radiates 
in  the  soft  X-ray  regime.  We  shall  see  that  the  accretion  disk  model 
motivates  us  to  identify  super  luminous  X-ray  sources  as  black  hole  accretion 
disks. 


1.2  Steady  State  Accretion  Disks 


Accretion  disks  account  for  the  black  hole’s  remarkable  efficiency  in  the 
conversion  of  gravitational  binding  energy  into  X-ray  luminous  energy.  As 
shown  in  section  1.1,  while  spherical  accretion  can  give  rise  to  conversion 
efficiencies  on  the  order  of  5%  of  the  rest  mass,  a  disk  about  a  rotating  black 
hole  can  have  conversion  efficiencies  of  up  to  40%  me2. 

In  general,  matter  flowing  toward  a  black  hole  is  expected  to  form  an 
accretion  disk  in  a  binary  system.  The  two  primary  mechanisms  by  which 
mass  can  be  injected  into  the  black  hole's  gravitational  well  are  Roche  lobe 
overflow  and  stellar  wind  capture.  Both  mechanisms  inject  mass  into  the 
black  hole's  gravitational  well  with  an  angular  momentum  per  unit  mass 
which  is  proportional  to  that  of  the  binary  system.  Some  of  the  injected  mass 
is  thrown  off,  and  some  is  accreted.  Since  the  accreted  mass  has  angular 
momentum  with  respect  to  the  black  hole,  the  mass  cannot  enter  the  hole 
directly,  but  must  form  a  disk  orbiting  in  the  plane  of  the  system.  Through 
viscous  stresses  the  matter  in  he  disk  looses  angular  momentum  and  spirals 
inward  until  it  looses  sufficient  angular  momentum  to  bring  it  within  the  last 
stable  orbit  predicted  by  general  relativity.  Only  then  can  it  enter  the  hole. 

The  goal  of  this  chapter  is  to  review  the  theory  of  accretion  disks  about 
black  holes  and  to  show  why  one  expects  an  accretion  disk  to  radiate  in  the 
soft  X-ray  regime  for  it  was  through  its  X-ray  radiation  that  A0620  was 
identified  as  a  compact  object. 


1 .2  Steady  State  Accretion  Disks 
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1.2. a  Structure  of  the  "Standard  Disk" 

In  can  be  shown  (Shakura  and  Sunyaev  (1973),  and  ST)  that  by  employing 
conservation  laws,  an  equation  of  state,  a  viscosity  law,  and  a  radiative 
transport  law  that  one  can  derive  the  basic  structures  and  spectra  or  an 
accretion  disk.  The  structures  and  spectra  of  the  disk,  in  turn,  are  functions 
of  three  parameters:  the  radius,  r,  the  mass  of  the  central  object,  M,  and  and 
the  mass  accretion  rate,  M.  For  fixed  values  of  these  three  parameters  the 
disk  can  be  conveniently  divided  into  three  regions. 

1.  An  inner  region  where  radiation  pressure  and  scattering  dominate. 

2.  A  middle  region  where  the  dominant  source  of  opacity  is  scattering, 
but  gas  pressure  dominates  radiation  pressure.  The  transition 
between  the  inner  and  middle  regions  occurs  at  the  radius  rmi  where 
Pg  =  Prad- 

3.  An  outer  region,  at  large  r,  where  gas  pressure  dominates  radiation 
pressure  and  the  opacity  is  dominated  by  free-free  absorption.  The 
transition  from  the  outer  to  the  middle  region  occurs  at  a 

radius,  rom»  where  Kff  ~  KeS. 
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1.2.b  Spectrum  of  the  "Standard"  Disk 


The  standard  model  predicts  disk  temperatures  on  of  the  order  T  *  107K. 
When  attempting  to  interpret  the  spectrum  of  the  disk,  it  is  important  to 
realize  that  the  temperature  calculated  in  the  standard  model  is  the 
temperature  at  midplane  (z=0).  The  observable  temperature  however  is  the 
temperature  at  the  effective  emitting  surface,  where  the  photons  which 
escape  to  infinity  are  generated.  We  expect  the  surface  temperature,  Ts  to  be 
less  than  the  midplane  temperature  except  when  the  disk  is  optically  thin,  in 
which  case  photons  generated  at  the  center  escape.  For  optically  thick  disks 
in  regions  where  absorption  dominates  scattering,  the  emission  has  a 
blackbody  spectrum  with  characteristic  temperature  Ts,  which  is  the 
temperature  of  the  disk  at  an  optical  depth  tfl  ~  1  below  the  surface  where  the 

average  escaping  photon  is  generated.  Therefore  the  emission  spectrum  is 
that  of  a  blackbody  in  the  outer  regions,  r  >  rom,  where  absorption  dominates 
scattering.  It  can  be  shown  that  in  these  regions,  Ts  <  105K. 

In  the  middle  to  inner  regions,  scattering  dominates  absorption  and  the 
blackbody  approximation  no  longer  holds.  Following  emission,  the  photons 
undergo  many  scatterings  which  can  best  be  characterized  as  a  random  walk. 
This  scattering  generates  a  modified  blackbody  distribution.  It  can  be  shown 
that  the  effective  blackbody  temperature  is  approximately 


Teff  ~  Ts 


'KffM/8 


(1.2.b) 


where  the  Rosseland  mean  opacity  is  used  to  evaluate  the  parentheses.  Since 
it  is  assumed  that  Kff  «  Kes  it  stands  that  Ts  »  Teff.  Since  Ts  represents  the 
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mean  temperature  of  the  region  emitting  the  photons,  one  can  see  that  the 
effect  of  scattering  is  to  raise  the  mean  energy  of  the  emitted  photons  (~kTs) 
above  the  blackbody  value  at  thermodynamic  equilibrium.  It  can  be  shown 
that  in  these  regions  Ts  is  on  the  order  of  106  -  107K  and  therefore  their 
emission  is  "harder." 

In  the  innermost  regions,  Comptonization  becomes  important.  In  this 
process,  photons  are  scattered  by  relativistic  electrons.  These  regions  are 
even  hotter,  reaching  T  <10®  K.  Figure  1.2.b  shows  a  typical  power  spectra. 


Sketch  of  the  power  spectra  emitted  by  two  model  diski  around  black  holes,  as 
calculated  bv  Shakura  and  Sunyaev  (1973)  without  taking  account  of  relauvisuc  conecuons  or  capture 
of  radiation  by  the  hole.  In  both  models  the  hole  was  assumed  to  be  nonrotaung,  so  the  inner  edge  of 
•he  disk  was  at  r,  —  6M  Model  ta)  corresponds  to  a  -  I0~\  M  ”  Vf0,  M  —  IO'8.Vf0  yr “ 1 , 
»  L FjJd  =  I0)S  erg  s'1;  model  ( b )  corresponds  to  a  -  I0~:  to  1  (spectrum  insensitive  to  a), 
M  -  M0,  M  -  l0"4Vo  yr" L  »  1034  erg  s" The  portion  of  the  spectrum  marked  "blackbody"  is 
generated  primarily  by  the  outer,  cold  region  of  the  disk  where  electron  scattering  is  unimportant.  The 
portion  marked  "  modified"  is  generated  by  the  middle  and  inner  regions  where  electron  scattering  is 
the  dominant  source  of  opacity  The  temperature  of  the  exponential  tail  is  the  surface  temperature  of 
the  mncnnosi  region.  (From  Novikov  and  Thome  ( 1973)  ] 


Figure  1.2.b:  Typical  accretion  disk  spectrum  (ST). 


1.3  Non-Steady  State  Accretion  Disks 


The  steady  state  results  presented  in  chapter  1 .2  are  insufficient  to  explain 
the  transient  nature  of  X-ray  sources  such  as  A0620.  Attempts  to  model 
these  types  of  sources  fall  into  two  categories:  disk  instability  models  and 
mass  transfer  instability  models. 

Within  the  disk  instability  picture  there  are  two  schools  of  thought. 
Mineshige  and  Wheeler  (1989),  Mineshige,  Tuchman  and  Wheeler  (1990) 
and  Huang  and  Wheeler  (1989),  HW  hereafter,  claim  that  it  is  possible  to  set 
up  temperature  and  density  waves  in  the  disk  as  it  accretes.  They  assume  a 
constant  low  mass  transfer  rate,  a  two  temperature  viscosity  parameter,  and  a 
temperature  dependent  opacity  to  account  for  atomic  and  molecular 
absorption.  This  leads  to  a  double  valued  relationship  between  temperature 
and  surface  density  and  the  model  undergoes  a  subcritical  pitchfork 
bifurcation  into  two  stable  surface  densities,  a  lower  value  Scold  and  a  higher 
value  Shot-  When  the  interior  of  the  slowly  accreting  disk  surpasses  Shot,  a 
heating  wave  propagates  outward  and  causes  what  HW  call  a  "minor" 
outburst.  If  the  heating  wave  propagates  far  enough  out  it  can  trigger  an 
inward  propagating  cooling  wave  which  causes  the  disk  to  accrete  strongly  in 
order  to  lower  the  surface  density  to  Scold  in  the  inner  regions.  This 
corresponds  to  what  HW  call  a  "major"  burst  such  as  the  1975  burst  of 
A0620.  One  can  choose  the  free  parameters  and  boundary  conditions  in  this 
model  to  give  good  agreement  with  the  A0620  outburst. 

Other  disk  instability  models  (Canizzo  et  al.  1982;  Lin  and  Taam  1984) 
are  reminiscent  of  models  for  dwarf  novae  outbursts  in  which  a  low 
accretion  rate  causes  the  density  of  the  disk  to  slowly  build  up  until  the  disk 
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becomes  optically  thick.  At  this  point  the  disk  undergoes  convection  which 
rapidly  transports  matter  onto  the  compact  object  creating  a  nova. 

Mass  transfer  instability  theories  attribute  the  nova  to  the  companion. 
Hameury  et  al.  (1986)  claim  that  hard  X-rays  from  accretion  can  induce  an 
instability  in  the  photosphere  of  the  companion.  At  the  inner  Lagrangian 
point,  hard  X-rays  can  penetrate  just  below  the  photosphere,  heating  it.  This 
X-ray  heating  expands  the  atmosphere  at  LI  causing  an  envelope  of  extended 
stellar  atmosphere  to  form.  When  the  distance  between  the  extended 
atmosphere  and  LI  reaches  a  critical  value,  matter  rapidly  streams  across  the 
LI  point  and  onto  the  accretion  disk.  This  increases  the  hard  X-ray  flux 
from  accretion  which  increases  the  X-ray  heating  of  the  companion  which 
causes  more  matter  to  enter  the  envelope,  etc.  This  positive  feedback 
continues  until  the  disk  becomes  optically  thick,  shielding  the  LI  point  and 
quenching  the  instability.  The  nova  then  decays  until  the  disk  exhausts  most 
of  its  mass.  The  authors  point  out  that  the  key  to  testing  this  theory  is 
observation  of  hard  X-rays.  Soft  X-rays  will  be  absorbed  by  the  photosphere 
and  merely  expand  the  corona,  not  the  atmosphere. 

All  of  the  above  models  predict  recurrent  novae  and  there  are  sufficient 
free  parameters  to  match  the  tentative  ~  70  year  period  between  bursts  for 
A0620. 


1.4  Candidates 

There  are  three  binary  x-ray  sources  for  which  there  is  strong  evidence  of 
their  being  black  holes.  In  each  case  the  mass  predicted  by  the  mass  function 
is  large  (M2  ~  3Mq),  there  is  flux  variability  on  short  time  scales  implying 
compact  size,  and  the  X-ray  luminosity  is  huge  (L  ~  1037  erg/s  -KHLq).  The 
three  candidates  are  CYG  X-l,  LMC  X-3,  and  A0620.  In  this  section  all  data 
such  as  mass  functions,  luminosities,  distances,  periods,  etc  are  taken  from 
the  concise  summary  given  by  McClintock  (2333)  unless  otherwise  stated. 

1.4. a  Cygnus  X-l 

Perhaps  the  best  known  of  all  black  hole  candidates,  CYG  X-l  was  the 
first  to  be  proposed  to  be  a  black  hole  in  a  binary  system.  The  compact  object 
has  a  mass  function  (equation  l.l.a.l)  of  .25  ±  .01  Mq,  which  taken  alone  is 
not  convincing  evidence  for  it  to  be  a  black  hole  binary. 

The  optical  companion  is  an  09.7  supergiant,  HDE226868.  If  we  assign  a 
typical  mass  to  HDE226868  of  33Mq  then  the  mass  of  the  compact  object 
cannot  be  less  than  7Mq.  However,  it  is  possible  that  a  companion  in  a  close 
binary  system  may  be  severely  undermassive  for  its  spectral  type 
(McClintock  2333;  Paczynski  1973;  van  den  Heuvel  1983)  and  thus  the  mass 
of  the  compact  could  be  substantially  less. 

A  second  way  of  determining  the  mass  comes  from  stellar  structure 
calculations.  While  the  supergiant  may  be  undermassive,  it  can  be  shown 
(ST,  Kippenhahn  1969)  that  the  mass  must  be  at  least  8.5  Mq  to  have  the 
correct  luminosity  at  a  distance  d  ~  2  kpc.  The  distance,  d,  can  be  found  by 
measuring  several  stars  within  about  a  degree  of  HDE226868  and  averaging 
the  absorption  and  then  using  this  absorption  to  adjust  the  reddening-distance 
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measurements.  Two  groups  have  done  this  and  concluded  that  "HDE  226868 
is  not  closer  than  1  kpc,  and  that  its  minimum  distance  is  about  2.5  kpc" 
(McClintock  2333;  Bregman  et  al.  1973).  Using  8.5  Mq  for  the  mass  of  the 
optical  star,  from  the  mass  function  we  obtain  a  mass  of  M2>  3.3  Mq. 

According  to  Paczynski  (1974),  more  conclusive  evidence  for  CYG  X-l 
to  be  a  black  hole  comes  from  the  lack  of  X-ray  eclipsing.  In  this  event, 
limitations  can  be  placed  on  the  inclination  of  the  orbital  plane  with  respect  to 
terrestrial  observers  and  nothing  has  to  be  assumed  about  the  mass  of  the 
supergiant.  Following  ST,  if  there  is  no  eclipsing,  it  must  be  true  that  (see 
figure  1.4.a)  cos  i  >  R/a 


Figure  1.4.a  No  eclipse  condition 


It  can  be  shown  that  the  no  eclipse  condition  translates  into  requiring  that 
M2sin  i  cos2  i  >  fR2/(aisin  i)2.  For  a  lower  bound  on  M2  one  substitutes  in 
(sin  i  cos2  i)max  =  2/3^3 . 
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One  can  relate  R2  to  its  luminosity  and  effective  temperature  to  get  an  overall 
result  of 


3V3  f  Lq 
“  8^aT4(aisin  i)2 


1()4+0.4(4.72-mv-BC+Av) 


(1.4. 


where  mv  is  the  apparent  magnitude  in  the  V,  BC  is  to  bolometric  correction 
and  Av  is  the  absorption  in  the  V.  Evaluating  this  expression  with  the 
parameters  of  CYG  X-l,  one  finds  that 

M2  >  3.4  M0  <2l^)2  (1.4.a.2) 

If  d  >  2  kpc  then  this  is  very  probably  a  black  hole;  if  the  distance  is  less,  then 
we  cannot  say  for  sure.  For  example,  if  d  =  1.5  kpc  then  M2=  1.9  Mq. 

If  the  supergiant  is  required  not  to  overflow  its  Roche  lobe  then  from  ST 

M2  >  6.5  Mq  (d  =  2  kpc).  From  all  of  these  considerations  one  is  forced  to 

conclude  that  the  mass  of  CYG  X-l  probably  exceeds  3  Mq  and  therefore, 

since  it  cannot  be  a  neutron  star,  it  is  probably  a  black  hole. 

1.4.b  LMC  X-3 

LMC  X-3  is  a  better  candidate  for  two  reasons.  It  has  a  mass  function 
f=2.3  ±  0.3  Mq  which  exceeds  the  maximum  neutron  star  mass  in  some 
models.  Secondly,  since  it  is  in  the  Large  Magellanic  Cloud,  its  distance  is 
known  to  be  d  =  55  kpc  with  small  uncertainty. 

Since  LMC  X-3  exhibits  no  x-ray  eclipse  one  can  use  formula  (1.4.a.l) 
with  the  appropriate  parameters  to  give  M2  >  6  Mq.  If  a  typical  mass  for  the 
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optical  companion's  spectral  type  is  attributed  to  it,  one  finds  M2  >  7  Mq 
(McClintock  2333;  Cowley  etal.  1983). 

1.4.C  A0620 

In  chapter  1.2  it  was  shown  that  an  accreting  compact  object  could  have  a 
huge  luminosity  in  the  soft  X-ray  region  of  the  spectrum.  It  is  difficult  to 
find  another  credible  physical  process  which  could  produce  this  type  of 
spectrum  and  luminosity.  The  luminosity  of  A0620  in  soft  X-rays  was  ~  1038 
erg  s'1.  Therefore  it  is  certainly  a  compact  object. 

The  strongest  evidence  that  A0620  is  a  black  hole  and  not  a  neutron  star  is, 
first  and  foremost,  its  optical  mass  function  which  is  3.18  ±  0.16  Mq  (MR). 
Thus,  in  the  unphysical  limit  of  a  zero  mass  optical  companion,  this  is  the 
minimum  mass  of  the  compact  object.  This  mass  is  much  greater  than  the 
maximum  mass  of  a  neutron  star  in  almost  all  current  theories.  From  the 
observed  mass  function  one  can  derive  an  absolute  3a  lower  limit  on  the  mass 
of  the  compact  object  of  M2  =  2.7  Mq. 

One  can  get  a  firm  lower  limit  by  assuming  that  the  optical  companion 
star  is  a  zero  age  main  sequence  star.  An  average  mass  for  a  K  dwarf  is  ~  0.5 
to  0.8  Mq.  Paczynski(1973)  and  van  den  Heuvel(1983)  have  shown  that 
binary  accreting  systems  may  have  undermassive  optical  companions. 
Therefore  a  minimal  mass  of  0.25  Mq  is  chosen  for  the  companion.  Since 
there  is  no  X-ray  eclipsing  a  credible  maximum,  i  =  85°,  is  chosen.  Choosing 
a  velocity,  vj,  of  the  K  star  which  is  3a  less  than  the  best  fit  value  one  obtains 
a  3a  lower  limit  of  3.20  Mq. 

One  can  find  a  realistic  mass  by  assuming  a  K5V  star  with  a  mass  of 
0.7Mq.  This  is  the  preferred  spectral  type  and  mass  for  the  K  dwarf .  In  this 
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case  the  3a  limit  is  M2  >  4.0  Mq.  If  one  further  requires  the  dwarf  to  fill  its 
Roche  lobe  and  the  value  suggested  by  the  amplitude  of  ellipsoidal  variations 
of  i  <  50°  is  used  for  the  inclination,  one  finds  M2  >  7.3  Mq  (3a). 

Another  method  of  placing  a  limit  on  the  mass  of  the  compact  object  has 
been  suggested  by  Johnston  et  al.  (1989).  They  observed  the  Doppler  shifts 
of  the  H«  and  Hg  lines  in  the  disk  and  fitted  them  to  an  accretion  model  and 

then  averaged  over  several  measurements  to  obtain  a  profile  similar  to  1.4.c. 


Velocity  (kms  ' ) 


Figure  1.4.c  Summed  profile  of  all  Ha  spectra  (Johnston  et  al.) 

Using  this  profile,  they  were  able  to  make  an  estimate  of  the  velocity  at  the 
outer  edge  of  the  disk.  Assuming  the  accretion  disk  is  in  the  binary  plane  and 
that  the  velocity  profile  is  given  by  circular  Keplerian  orbits,  v(r)=VGMx/r, 
the  outer  edge  represents  the  minimum  velocity.  If  one  requires  that  the 
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outer  edge  of  the  disk  be  within  the  Roche  lobe  of  the  black  hole  then,  solving 
numerically,  the  absolute  3a  lower  limit  is  M2  >  4  Mq  if  Mi  =0.25  Mq.  If 

the  mass  of  the  companion  is  0.7  Mq,  then  the  absolute  3a  lower  limit  is  M2  ^ 

9.3  Mq. 

In  summation,  we  see  that  there  are  several  strong  black  hole  candidates. 
A0620  is  the  most  convincing  of  the  three. 


Chapter  2:  Data  Reduction 


In  this  chapter,  I  discuss  the  various  computer  codes  which  were  used  to 
reduce  the  raw  data  from  the  SAS-3  satellite.  I  also  explain  how  errors  were 
assessed  on  the  individual  measurements.  I  examined  data  from  observations 
of  the  Crab  Nebula  to  obtain  objective  criteria  for  the  selection  and  rejection 
of  fitted  data  points  and  obtain  elevation  dependent  response  functions  for  the 
slat  collimators.  Finally,  I  discuss  the  probable  systematic  error  induced  into 
the  data  set. 

2.1  Fitting  the  Data 

The  software  packages  used  to  reduce  the  data  are  nbin ,  aspect,  clean  and 
fit  all  of  which  were  written  by  Richard  Brezina.  The  program  nbin  decodes 
the  data  from  the  SAS-3  data  archives  and  places  it  in  tabular  format.  Aspect 
is  used  to  determine  satellite  position,  field  of  view  for  the  detector  of 
interest  and,  where  applicable,  the  transmission  function  for  the  source  at  a 
given  moment  in  time.  Clean  separates  out  source  scans  in  the  data  stream 
output  by  aspect  and  attaches  a  portion  of  the  data  stream  before  and  after  the 
scan  for  background  calculations.  It  also  insures  that  there  were  occultations 
neither  by  the  Earth  nor  the  Moon  and  checks  for  the  South  Atlantic 
Anomaly,  any  of  which  is  grounds  to  discard  the  scan.  Fit  then  generates  the 
relevant  statistics.  For  a  given  orbit  each  routine  is  performed.  I  will  only 
discuss  fit  in  detail. 

The  fit  program  fits  scanning-mode  data  to  the  triangular  transmission 
function  generated  by  aspect  as  shown  in  figure  2.1  (next  page).  It  tries  to 
correct  for  satellite  timing  errors  by  matching  up  the  maxima  of  the  least 
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squares  fit  of  the  aspect  calculation  and  the  actual  scan.  The  fit  is  performed 
using  the  curfit  routine  given  in  Bevington  (1969).  The  routine  returns  a 
fitted  intensity,  Io,  a  fitted  count  rate  for  each  measurement  time,  and  the  y} 
for  each  fit  with  typical  reduced  x2  =  1.0.  We  define  the  error  in  Io  by 


V2,counts+(scan  duration)x(average  background  rate)  n 
Cstat  -  ,  (“• 

J  (transmission  function)dt 


where  ^counts  is  the  total  counts  in  the  scan  and  transmission  function  is  the 
percent  transmission  of  the  source  during  each  time  bin.  Thus  (2.1.1) 
represents  the  Poisson  error  associated  with  the  source,  but  not  errors  in  the 
transmission  function  determined  by  aspect.  In  chapter  2.2, 1  address  this 
problem 


I  terr  , 


figure  2.1:  Typical  data  scan.  Fit  adjusts  the  temporal  position 
of  the  aspect  curve  by  y}  minimization 


2.2  Systematic  Corrections 


In  an  effort  to  assess  the  systematic  errors  in  the  intensities  derived  from 
A0620  data,  we  analyzed  similar  data  on  the  Crab  Nebula.  Slice  the  Crab  is  a 
steady  source  on  all  time  scales  longer  than  ~  1  second,  any  systematic 
deviations  of  the  computed  intensity  away  from  a  constant  value  is  a  symptom 
of  error  which  may  also  affect  the  accuracy  of  the  computed  A0620 
intensities.  Using  the  Crab  as  our  standard  "candle"  we  derived  objective 
empirical  criteria  for  selection  and  correction  of  the  A0620  intensities.  We 
also  derived  a  response  function  for  each  channel  which  corrects  for  the 
systematic  variation  of  intensity  with  elevation  angle.  In  the  following 
section,  when  we  speak  of  data  points  we  are  refering  to  the  fitted  intensity 
generated  by  fit  using  the  raw  SAS-3  count  data.  Using  the  center  slat  A 
channel  as  an  example,  we  describe  how  the  tasks  were  accomplished  in  seven 
steps. 

1.  Purification 

We  discard  any  data  points  where  A0620  and  Crab  are  simultaneously  in 
the  channel's  field  of  view.  This  can  be  done  trivially  by  comparing  crossing 
times  and  sweep  durations.  Figure  2.2a  illustrates  the  raw  Crab  data  and 
figure  2.2b  shows  the  data  after  purification. 

2.  Remove  Poor  Fits  and  Excessive  Timing  Errors 

We  discard  points  for  which  the  reduced  x2  °f  the  triangle  fit  to  the 
scanning-mode  data  is  greater  than  two.  We  do  this  since  we  are  ultimately 
going  to  try  to  derive  an  elevation  dependent  response  function  based  on  the 
data  and  we  only  want  to  use  our  best  data. 


Center  Slat  A  Channel  Crab  Counts/Sec 
100  200  0  100  200  0  100  200 
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Figure  2.2:  (a)  The  unprocessed  intensity  fits. 

(b)  The  data  after  removing  scans  when  A 0620  was  also  in  the  field  of  view 

(c)  Points  with  y}  >  2  and  timing  errors  >  ±8.64  seconds  shown  as  stars. 
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We  also  remove  points  with  excessive  timing  errors.  By  excessive  we 
mean  greater  than  ±  8.64  seconds  =  1CM  day.  The  reason  behind  this 
discarding  is  that  large  timing  errors  represent  data  points  with  questionable 
spacecraft  attitude  information.  In  figure  2.2c,  points  removed  by  these 
screening  criteria  are  shown  as  stars. 

3.  Removal  of  Other  Less  than  Desirable  Points 

All  of  the  reduction  software  is  designed  to  be  used  on  SAS  counting  rates 
obtained  in  the  scanning  mode.  A  typical  scan  over  a  source  lasts  for  15  -  45 
seconds.  Counts  obtained  in  the  pointing  mode  are  not  handled  at  all  by  these 
routines.  Ultra-fast  sweeps  (At  <  8  sec)  are  not  handled  well.  Therefore  we 
remove  all  points  derived  from  scans  with  duration  not  in  the  range  8  <  At  < 
50  seconds. 

Another  source  of  error  is  high  background  count  rates.  For  background 
determinations  we  use  data  from  100  seconds  before  and  after  each  scan.  If 
another  smaller  source  enters  the  detector  field  of  view  during  these 
background  margins,  it  will  perturb  the  apparent  background  rate.  Since  we 
subtract  the  background  count  rate  from  the  source  count  rate,  this  can  lead 
to  inaccurate  results.  We  therefore  define  the  signal  to  noise  ratio  (SNR)  to 
be  Io/b  where  b  is  the  background  count  rate  and  Io  is  the  fitted  intensity.  For 
the  high  count  rate  A  and  B  channels  we  allow  SNR  >  5.  For  the  C  channel, 
which  has  a  lower  count  rate,  we  are  forced  to  accept  SNR  >  2.  For  the 

extremely  low  count  rate  D  channel  we  accept  any  point  with  SNR  >  1 .4. 

Finally,  we  discard  data  for  which  a/Io  >  Such  a  point  is  of  little 

statistical  significance  since  its  error  bars  are  so  large.  Figure  2. 2d  shows 
these  other  less  than  desirable  points  as  stars. 
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4.  Remove  Points  Outside  the  FWHM 

As  can  be  seen  from  figure  2.2d,  points  outside  ±35°  are  unreliable.  We 
therefore  remove  these  points. 

5.  Window  Elevation  Angle  Corrections 

The  beryllium  windows  on  the  slat  detectors  absorb  soft  X-rays.  If  a 
photon  crosses  the  window  at  an  angle,  0,  the  effective  thickness  of  the 
window  is  increased  by  the  factor  sec9  which  causes  increased  absorption. 
To  correct  for  this  effect  we  calculate  the  window  transmission  at  a 
particular  angle,  tx(0).  We  then  relate  this  value  to  the  transmission  to  the 
value  at  0  =  0  by  the  percent  transmission  formula 

Emax 

tx(0)  =  ^o)  Jexp[-TjiBesec0]exp[-o(E)NH]E-«dE  (2.2.1) 
Emin 

where  a(E)  is  the  energy  dependent  interstellar  photoelectric  absorption 
cross  section,  Nh  is  the  hydrogen  column  density,  a  is  the  spectral  index,  jJ.Be 
is  the  mass  absorption  coefficient  for  beryllium,  T  is  the  thickness  of  the 
window,  Emin  and  Emax  are  the  channel  boundaries  of  the  detector  in 
question  and  we  have  turned  the  transmission  at  an  angle  into  a  percentage  of 
the  transmission  at  0  =  0  by  dividing  by  tx(0),  the  window  transmission 
function  at  0°.  From  Woo(1991)  the  photon  distribution  in  the  energy 
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Figure  2.2:  (d)  Points  with  SNR  <  .5,  c/Io  >  -5  and  not  from  scans  depicted  as  stars. 

(e)  The  data  before  removing  points  outside  the  FWHM. 

(f)  The  data  after  window  correction. 
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bands  of  interest  can  be  fit  by  using  Nh  =  3.45x1021  particles/cm1 2  and  a  = 
2.1.  The  data  are  then  corrected  as 

io’  =  io/tx(e) 

o'=  |l0'  =  <J/tx(0)  (2.2.2) 

where  Io  is  the  fitted  intensity  and  a  is  the  Poisson  error  associated  with  it.  In 
equation  (2.2.2)  we  maintain  the  same  percentage  error  bars  i.e.  if  we  knew 
Io  to  within,  say,  10%  before  the  correction,  we  still  know  it  to  within  10% 
after  correction.  We  maintain  relative  error  bars  because  we  are  correcting 
the  data  for  systematic  effects.  Later  in  this  chapter,  I  will  discuss  assessing 
additional  systematic  errors  based  on  these  systematic  effects. 

Figure  2.2e  shows  the  data  after  the  window  correction.  One  can  see  that 
there  is  a  systematic  trend  with  respect  to  elevation.  We  now  proceed  to 
correct  for  this  effect. 


6.  Generate  the  Systematic  Correction  Function, 

We  now  generate  the  systematic  correction  function  due  to  elevation 
angle,  0  through  an  iterative  process: 

1.  Fit  the  data  to  a  3rd  order  polynomial  y  =  a  +  b@  +  c02  +  d03 4 5. 

2.  Subtract  the  quadratic  and  cubic  terms  c02  +  d03. 

3.  Set  avg  =  the  average  of  the  new  (nearly)  linear  data  set. 

4.  Refer  back  to  the  original  (pre-polynomial  fit)  data  set  and 
remove  any  points  which  are  more  than  3a  away  from  the 
average,  avg,  calculated  in  step  3. 

5.  Repeat  from  step  1  until  no  points  are  discarded  in  step  4. 
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Chap 

Figure  2.2f  shows  the  corrected  data  from  step  2.  The  reduced  x2  for  the 
linear  fit  is  1.085. 
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(g)  Corrected  Data 
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Figure  2.2:  (g)  The  corrected  intensity  fits  with  error  bars. 

(h)  The  error  dispersion.  The  curve  is  the  minimum  y}  fit  to  a  Gaussian. 
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We  now  construct  the  angular  response  function  which  we  normalize  by 
the  average  intensity  to  make  it  independent  of  source.  For  each  elevation 
angle  0,  we  associate  a  correction  factor  f(0)  given  by 

f(0)  =  a  +  b0  +  c02  +  d03 

Io"  =  lo'f(0) 

o"  =  o’f(0)  (2.2.3) 

where  Io'  and  o'  are  the  results  from  section  5  above.  We  now  apply  our  new 
response  function  to  the  pre-polynomial  fit  data.  The  reduced  x2  is  1 .08  (446 
d.o.f.)  and  the  new  data  is  shown  in  figure  2.2g. 

For  completeness,  the  following  table  summarizes  the  results  of  this 
section  (note:  all  coefficients  are  given  for  0  given  in  units  of  -k/2  <  0  <  ti/2. 

Table  1 


Channel  Crab  Correction  Coefficients  and  Averages 


Coefficient 

Channel 

a 

b 

c 

d 

avg 

Center  A 

1.000 

-0.084 

-0.328 

-0.615 

140.369756 

Center  B 

1.000 

-0.131 

0.131 

-0.378 

30.964078 

Center  C 

1.002 

-0.105 

0.346 

0.245 

20.613074 

Center  D 

2.165 

-0.769 

13.180 

12.544 

.694907 

Left  A 

1.000 

0.004 

-0.383 

-0.297 

137.685564 

Left  B 

1.001 

0.122 

0.101 

-2.261 

35.134279 

Right  A 

1.000 

0.058 

-0.630 

-0.701 

141.724995 

Right  B 

1.003 

0.138 

0.813 

-1.089 

22.635886 
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It  should  be  noted  that  the  result  for  the  D  channel  is  highly  questionable 
due  to  the  low  source  count  rates  and  correspondingly  large  relative  Poisson 
error. 

7.  Additional  Systematic  Error 

The  question  of  whether  there  should  be  additional  error  added  to  the 
Poisson  error  is  a  delicate  one.  Figure  2.2.h  is  a  histogram  of  the  number  of 
center  slat  A  points  that  are  a  given  number  of  o  away  from  the  average 
given  in  the  table  above.  Superimposed  over  the  histogram  is  a  Gaussian  fit. 
The  fit  has  a  reduced  y}  of  0.898  which  translates  into  a  confidence  of  68% 
and  so  is  a  fairly  good  fit.  If  the  only  source  of  error  were  Poisson 
fluctuation  we  would  expect  an  approximate  Gaussian.  Since  the  distribution 
of  error  is  nearly  Gaussian  we  can  conclude  that  the  residual  systematic 
errors  are  as  small  as,  or  smaller  than,  the  Poisson  error  and  therefore,  to 
first  order,  negligible.  The  resulting  data  set  seems  to  be  as  unbiased  and 
Poissonian  as  possible  and  thus  no  further  corrections  are  made.  The  net 
results  from  this  analysis  are: 

1.  We  have  determined  a  set  of  objective  criteria  through  which 
to  screen  the  data  set. 

2.  We  have  written  an  elevation  dependent  response  function 
based  on  real  data  which  corrects  for  systematic  defects  in 
each  of  the  channels. 

3.  We  now  have  confidence  in  the  corrected  data  set.  Any 
morphology  or  periodicity  observed  in  the  data  is  likely  to  be 
true  and  not  the  result  of  some  sort  of  processing  artifact  or 
systematic  defect. 


Chapter  3:  Application  to  A0620 


The  attribution  of  errors  to  the  data  is  the  least  certain  part  of  the 
reduction  scheme.  There  are  systematic  errors  which  may  be  as  large  or 
even  larger  than  the  Poisson  errors  associated  with  the  count  rate.  In 
particular,  the  systematic  detector  defects  due  to  elevation  angle  are  a  major 
subject  to  address.  In  chapter  2, 1  showed  how  the  Crab  Nebula  was  used  to 
evaluate  these  systematic  defects  and  correct  for  them.  In  this  chapter  I  use 
these  results  in  the  analysis  of  the  A0620  data. 

A  paramount  concern  is  to  retain  as  much  data  as  possible  without 
compromising  our  confidence  in  it.  In  this  spirit,  the  following  eight-step 
process  was  used.  Figure  3.2  illustrates  each  step  as  carried  out  on  the  center 
slat  A  detector  (1.3-5  keV).  Figure  3.3  shows  the  final  results  for  the  A  and 
B  detectors. 

L  Exiracias  mu£h-data  as  possible  from,  rise* 

One  of  the  most  interesting  sections  of  data  was  that  obtained  during  the 
initial  rapid  rise  in  the  intensity  from  background  to  peak  intensity 
immediately  after  the  discovery  of  A0620.  Unfortunately,  the  position  of 
A0620  was  initially  not  known  to  a  high  degree  of  accuracy.  In  the  rush  to 
locate  and  identify  the  source,  many  satellite  aspect  changes  were  made  and  as 
a  result,  it  is  difficult  to  derive  accurate  position  data.  There  are  two  major 
errors:  timing  errors,  and  changes  in  the  spin  rate  during  a  source  scan. 

A  timing  error  occurs  if  the  source  is  not  in  the  field  of  view  when  the 
aspect  program  predicts  that  it  will  be.  Small  timing  errors,  on  the  order  of 
a  second,  are  easily  corrected  (see  chapter  2.1).  Large  timing  errors,  on  the 
order  of  tens  of  seconds,  are  more  difficult  to  correct.  There  may  be  several 
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sources  in  the  scan,  but  we  are  required  to  pick  the  correct  source.  This  is 
difficult  late  in  the  decay  phase  of  the  nova,  however  in  the  early  part  of  the 
light  curve,  especially  during  the  rise,  A0620  is  much  brighter  than  anything 
else  in  the  X-ray  sky  (on  the  order  of  500  Crabs!)  We  can  therefore  slide  the 
aspect  triangle  response  function  "by  hand"  into  the  correct  time  interval  and 
then  ca"v  out  the  fit  described  in  2.1  while  maintaining  a  high  degree  of 
confidence  in  the  data. 

Dealing  with  spin  rates  that  are  changing  during  the  source  scan  is  not  so 
simple.  First,  the  light  curve  is  no  longer  a  triangle  as  shown  in  figure  3.1. 
Second,  frequently  the  spin  information  is  not  contained  in  the  orbital 
information  data  files  and  aspect  uses  the  wrong  spin  rate.  This  results  in 
triangle  response  functions  with  the  wrong  base  width  (fig  3.1). 


counts  peak 


wrong  base  time 


Figure  3.1:  Spin  rate  change  during  a  source  scan.  The  source  scan  is  much  wider  than 
the  aspect  prediction.  Also  the  shape  is  not  triangular. 
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In  an  attempt  to  recover  data  afflicted  by  these  problems,  we  assume  that 
the  elevation  information  is  correct  and  use  the  peak  of  the  aspect  curve  with 
the  peak  of  the  counts  to  obtain  the  intensity  according  to  the  formula 

peak  counts  -  average  background  . 

peak  transmission  function  ^  *  ' 


where  the  average  background  is  obtained  by  examining  the  data  on  either 
side  of  the  scan.  We  associate  Poisson  errors  with  each  of  these  points  in  the 
usual  fashion 

Vpeak  counts  +  avg  background 
a  “  peak  transmission  function  ^  '  ' 


Empirically,  these  points  agree  well  with  the  other  data  points  generated  by 
the  fitting  so  we  maintain  a  high  degree  of  confidence  in  the  data  quality. 
Figure  3.2a  shows  raw  data  and  3.2b  shows  the  same  data  with  a  corrected 
rise  time. 

2.  Eliminating  Data  Contaminated  bv  X-ravs  from  the  Crab  Nebula 

The  slat  detectors  have  a  very  large  field  of  view  (2°x70°  FWHM  for  the 
center,  l°x70°  for  the  left  and  right  slats).  Other  sources  may  be  in  the  field 
of  view  at  the  same  time  as  A0620  and  affect  the  intensity  fitting.  The  only 
source  of  consequence  near  A0620  (RA  95.0°,  DEC  0.0°)  is  the  Crab  Nebula 
(RA  82.9°,  DEC  22.0°).  It  so  happens  that  on  rare  occasions  both  sources  are 
in  the  field  of  view  at  the  same  time.  When  this  happens,  we  find  the 
percentage  overlap  of  the  triangle  response  functions  and  multiply  this  by  the 
average  crab  intensity  for  that  detector  (average  intensity 
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Figure  3.2:  (a)  Unprocessed  intensity  fits. 

(b)  Data  after  improvement  of  the  initial  rise. 

(c)  After  eliminating  data  contaminated  by  the  Crab  Nebula. 
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determined  in  Chapter  2.2).  We  then  subtract  these  counts  from  the 
measured  counts  of  A0620  and  these  are  our  new  counts.  Intensities  and 
Poisson  errors  are  then  calculated  in  the  usual  fashion  using  these  new  counts. 
Again,  the  corrected  data  tracks  well  with  the  untainted  data  and  so  we  retain 
a  high  degree  of  confidence  in  our  correction  algorithm.  The  purified  data 
are  shown  in  3.2c. 

3.  Full  Width  at  Half  Maximum 

According  to  the  description  of  the  instrumentation  on  SAS-3  (Taylor, 
1975)  the  FWHM  of  the  slat  detectors  is  ±  35°.  From  the  observations  of  the 
Crab  Nebula  (Chapter  2),  I  find  that  outside  this  range  the  data  cannot  be 
corrected  reliably.  Therefore  I  remove  all  data  outside  the  FWHM.  Figure 
3. 2d  shows  the  corrected  data  for  the  center  slat  A  channel. 

4.  Remove  Timing  Errors 

Next  we  remove  data  with  timing  discrepancies  greater  that  ±8.64  seconds 
(10*4  days).  In  analyzing  the  initial  rapid  rise  we  expended  a  great  deal  of 
effort  to  use  even  this  type  of  data  However,  as  the  nova's  intensity 
exponentially  decays,  it  becomes  harder  and  harder  to  discern  A0620  from 
other  sources.  Fortunately,  there  are  not  many  scans  that  have  excessive 
timing  errors  and  we  do  not  have  to  discard  very  many  data  points.  Figure 
3.2e  shows  discarded  data  as  stars. 

5.  Other  Problem  Data 

Data  are  discarded  for  several  other  reasons  which  were  explained  in 
more  detail  in  Chapter  2.  The  following  is  a  brief  summary. 

If  the  data  come  from  a  pointing  mode  observation,  the  scanning  mode 
analysis  package  is  no  longer  appropriate  and  the  data  are  discarded.  If  the 
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Figure  3.2:  (d)  The  data  after  removing  points  obtained  outside  the  FWHM. 

(e)  Points  with  timing  errors  shown  as  stars. 

(f)  Points  with  SNR  <  .5,  c/Io  >  .5  and  not  from  scans  depicted  as  stars. 
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data  come  from  a  scan  with  a  very  high  spin  rate,  we  find  (from  the  Crab 
data)  that  it  is  also  poorly  fit,  and  must  be  discarded. 

We  also  require  that  the  signal  to  noise  ratio  be  within  acceptable  limits 
(empirically  2  for  A  and  B  channels  which  have  high  count  rates).  For  the  C 
and  D  channels  have  such  low  count  rates  that  this  test  is  not  appropriate. 
Finally,  we  require  that  the  data  be  meaningful  in  that  a/Io  >  -5.  We  only 
require  this  of  the  A  and  B  channel  data,  the  C  and  D  channels  being  so  low 
that  this  criterion  loses  relevance.  Figure  3.2f  shows  the  data  removed  by 
these  tests  as  stars. 

6.  Correct  for  Window  Transmission  Angle 

As  explained  in  Chapter  2,  the  beryllium  window  is  a  source  of  low 
energy  X-ray  absorption.  Increased  elevation  angles,  0,  correspond  to 
effective  window  thicknesses  increased  by  a  factor  of  secB.  From 
Woo(1991)  we  make  the  window  correction  using  an  average  hydrogen 
column  density  of  Nh  =  3.5xl022  particles/cm2  and  an  average  spectral 
index  of  a  =  6.43.  Figure  3.2g  depicts  the  center  slat  A  channel  data  after 
window  correction. 

7.  Correct  Systematic  Errors  Due  to  Elevation 

In  Chapter  2  I  derived  a  systematic  correction  factor  due  to  elevation  for 
each  detector.  For  each  elevation  angle,  0,  there  is  associated  a  percentage 
correction  f(0).  Figure  3.2h  depicts  center  slat  channel  A  after  this 
correction  is  carried  out. 
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Figure  3.2:  (g)  Data  after  window  correction. 

(h)  Data  after  elevation  angle  dependent  response  functions. 

(i)  The  final  data  set.  Note  improvement  over  fig.2.2a. 
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8,  Remove  Excessive  Y2 

Lastly  we  remove  data  points  for  which  the  reduced  x2  >  3.0.  This 
removes  points  which  were  poorly  fit  by  the  fitting  routine.  We  now  have  a 
high  degree  of  confidence  in  the  data  set  and  any  deviations  from  the  norm 
are  probably  significant  features  of  the  data  and  not  artifacts  of  the 
processing.  The  final  light  curve  is  shown  in  figure  3.2i  and  is  quite  an 
improvement  on  3.2a. 
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Figure  3.3:  (a)  The  Center  Slat  A  channel  final  data  set 


150  200  250  300 

SAS  Day  -  JD  -  2442500 


lication  to  A0620 


41 


OOC  002  001  0 

oas/s^unoo  a  tauu^no  0290V 

Figure  3.3:  (b)  The  Center  Slat  B  channel  final  data. 


140  160  180  200  220  240  260  280  300 

SAS  Day  =  JD  -  2442500 


42 


Chapter  3:  Application  to  A0620 
oas/siunoo  lautnaqo  g  0290V 
OOS  002  00T  0 


0009  OOOt  0002  0 

09S/s^unoo  lauueqo  y  0290V 


Figure  3.3:  (c)  Center  Slat  A  and  B  channel  plotted  together. 
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After  reducing  the  data,  and  attributing  the  proper  error  to  each  datum, 
we  can  proceed  to  extract  the  desired  information  about  the  presence  of 
absence  of  periodicities.  In  particular,  we  are  interested  in  whether  there  is 
evidence  for  the  observed  7.75h  optical  period  found  by  MR.  We  also  seek 
to  confirm  the  7.8d  X-ray  period  found  by  MB.  The  7. 8d  period  is  of 
particular  interest  since  MB  used  the  same  data,  but  without  the  improved 
data  reduction.  We  will  show  that  the  improved  reduction  yields  a  data  set  so 
"clean"  that  it  allows  us  to  place  strict  upper  limits  on  the  amplitude  of 
oscillations 

We  examine  the  long  term  scanning  mode  data  from  the  slat  collimator 
detector  by  the  method  of  phase  dispersion  minimization. 

4.1  Phase  Dispersion  Minimization 


The  idea  of  phase  dispersion  minimization  is  to  fold  the  data  at  all  possible 
periods,  11,  within  the  range  Flmin  to  Flmax-  Following  Stellingwerf(1978) 
we  define  the  variance  to  be 


a2  = 


I(*i  -  x  )2 
N-l 


(4.1.1) 


where  the  sum  is  from  1  to  N,  xi  is  the  i^1  datum  and  we  have  a  total  of  N  data 
points.  For  any  subset,  or  bin,  of  {xj}  we  define  the  sample  variance  s2 
exactly  as  in  (4.1.1)  except  that  N  is  replaced  by  n,  the  number  of  elements  in 
the  subset. 
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Suppose  that  we  have  M  subsets  such  that  the  variance  of  the  j*  subset  is  sj 

and  the  number  of  elements  in  it  is  nj.  We  define  the  overall  variance  for  all 
of  the  subsets  to  be 


S(nj  -  l)Sj2 

s2  =  Inj-M 


(4.1.2) 


Suppose  that  we  choose  these  subsets  to  be  phase  bins.  Choose  a  period  n 
and  assign  each  xi  to  a  phase  0i  =  ti  moddl)  where  ti  is  the  time  of 
measurement  xi.  Then  we  can  divide  the  set  { x }  into  M  bins  where  all  the 
members  of  bin  j  have  similar  <j>i.  Define  the  statistic 


0  =  s2/o2.  (4.1.3) 

If  we  have  chosen  the  correct  period,  then  s2  will  be  very  small  compared  to 
neighboring  periods  and  0  will  be  a  local  minimum.  If  FI  is  not  the  true 
period  then  s2  *  o2  and  0=1.  We  therefore  seek  to  minimize  0. 

In  practice  we  pick  Nc  covers  on  Nb  bins  for  a  total  of  M  =  NcNb  subsets. 
We  then  add  these  subsets  into  Nb  bins  of  width  Nc. 

To  analyze  the  statistical  certitude  of  any  observed  periodicity  we  follow 
Stellingwerf(1978)  and  note  that  ©  is  the  ratio  of  the  variances  of  two  subsets 
of  the  data  set,  that  of  the  total  set  and  that  of  the  bins,  and  therefore  obeys  an 
F  distribution  with  £nj  -  M  and  N  - 1  degrees  of  freedom.  It  is  convenient  to 
have  F  >  1  so  define  F  =  0*1  when  0  <  1  and  F  =  ©  otherwise.  For 
convenience,  I  will  assume  that  0<l(if0>l  replace  0  by  0*1). 
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We  can  apply  the  two-sided  test  to  determine  the  probability,  P,  that  a 
given  value  of  0  is  due  to  random  fluctuations.  This  "significance"  is  twice 
the  area  of  the  F  distribution  above  0*1.  The  significance  P  is  given  by 

F(p/2,N-1,  Xnj  -  M)  =  1/®  (4.1.4) 


we  then  find  P  by  consulting  an  F  table. 

Stellingwerf  points  out  that  this  can  be  simplified  by  noting  that  if  N  is 
large  (>100)  then  the  variance  of  our  sample  approaches  the  true  population 
variance.  In  that  case  we  can  use  the  y}  test 


X2(p/2,N-l)  =  (Inj  -  M)0 


(4.1.5) 


Now  P  is  twice  the  area  of  the  y}  distribution  below  (Xnj  -  M)0. 

The  above  discussion  places  no  upper  limit  on  periodicity  amplitude.  In 
chapter  4.2,  I  will  discuss  how  to  place  upper  limits  on  the  periodicity 
amplitude  using  the  observed  value  of  0. 

Since  the  data  are  rapidly  decaying  and  are  being  sampled  periodically,  it 
is  necessary  to  detrend  the  data.  Failure  to  do  so  can  cause  false  periods  to 
appear  and  also  cause  the  effects  of  sampling  to  grow  so  severe  as  to  mask  out 
any  periodicity  which  may  be  in  the  data  set  (Remillard,  1990).  To  do  this  I 
subtract  a  decaying  exponential  function  from  the  data  and  then  normalize 
according  to  the  formula 


yi  -  Iie-^ 

Iie-^ 


yi  = 


(4.1.6) 
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where  yi  is  the  data  point  at  time  t,  and  Ii  and  x  are  the  least  squares  fit 
parameters  for  the  data  set.  This  removes  the  long  term  trend  in  the  data  and 
allows  us  to  examine  departures  from  the  mean  value,  and  to  search  for 
periodicity. 

In  Chapter  4.2, 1  will  apply  these  techniques  to  the  A0620  data  set  in  order 
to  derive  limits  on  the  X-ray  periodicity. 


4.2  Results 


In  Chapter  4.1  we  discussed  the  development  of  methods  for  analyzing  the 
data.  In  this  chapter  we  apply  those  methods  to  data  from  the  A  and  B 
channels  of  the  slat  collimator  detectors.  We  chose  these  channels  because 
they  provided  the  largest  data  sets  (i.e.  the  most  points).  We  also  examine  the 
hardness  ratio  B/A.  In  section  4.2.1,  we  detrend  the  data  in  preparation  for 
the  analysis.  In  section  4.2.2,  we  describe  the  Monte  Carlo  simulation  we 
performed  to  analyze  the  data  and  present  our  results.  In  chapter  4.2.3,  we 
confirm  MB's  data  and  examine  their  7.8d  period  claim. 

4.2.1  Detrending  the  Data 

We  detrend  the  center  slat  A  channel  decay  data  as  follows.  We  ignore  the 
initial  burst.  We  fit  the  initial  decay  from  SAS  day  137  to  175  with  an 
exponential.  We  ignore  the  "hump"  between  SAS  days  175  and  200  and  then 
fit  another  exponential  to  the  decay  from  SAS  day  200  to  260.  Since  we  are 
looking  for  periodic  deviations  from  the  mean  and  not  singular  events,  we 
also  remove  the  hardness  "flare"  at  SAS  day  ~  215.  We  concatenate  the  two 
fits  together  for  periodicity  searches  and  we  name  this  section  of  the  data  the 
decay  phase.  Figures  4.2.1a-c  show  this  process  for  the  A  channel  where  the 
curves  represent  fits  and  the  error  bars  represent  the  data  and  uncertainty. 

Detrending  the  final  ramp  phase  of  the  data  is  difficult.  There  are  few 
natural  breaking  points  by  which  to  change  models  and  our  breaking  points 
are  admittedly  rather  arbitrary.  We  detrended  the  data  from  SAS  day  275  to 
285  by  fitting  a  line  to  it.  We  remove  the  average  of  the  data  between  days 
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Figure  4.2.1:  (a)  Channel  A  fit  for  SAS  days  137  -  175.  Curve  is  best  fit  exponential. 

(b)  Channel  A  fit  for  SAS  days  200  -  255.  Curve  is  best  fit  exponential 

(c)  Combined  percent  deviations  from  models. 
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Figure  4.2.2:  (a)  Channel  B  fit  for  SAS  days  137  -  175.  Curve  is  best  fit  exponential. 

(b)  Channel  B  fit  for  SAS  days  200  -  255.  Curve  is  best  fit  exponential 

(c)  Combined  percent  deviations  from  models. 


%  Deviation  from  Model  Hardness  Ratio  B/A 
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140  160  180  200  220 
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Figure  4.2.3:  (a)  Hardness  ratio  fit  for  SAS  days  135  -  240.  Curves  are  best  fits, 
(b)  Combined  percent  deviations  from  models. 
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285  and  300  and  between  301  and  307.  We  then  fit  the  remaining  data  from 
day  208  to  338  by  a  line.  Figure  4.2. Id  shows  the  fits  as  curves  and  the  data 
uncertainty  as  error  bars.  Figure  4.2.2e  shows  the  resultant  percent 
deviation  from  the  model.  In  the  analysis  section,  when  we  refer  to  the 
"ramp",  "MB"  and  "improved  MB"  sections  of  the  data  we  will  be  referring 
to  the  sections  from  day  308  to  338,  day  285  to  322,  and  day  275  to  308 
respectively. 

We  detrend  the  center  slat  B  channel  in  an  analogous  manner  to  the  center 
slat  A  channel  and  we  show  the  results  in  figure  4.2.2.  As  before,  the  curve 
represents  the  model  and  the  error  bars  represent  the  data  uncertainties.  We 
do  not  fit  the  B  channel  after  SAS  day  240  since  it  is  at  the  background  level 
and  is  no  longer  informative. 

We  detrend  the  hardness  ratio  B/A  by  fitting  with  an  exponential  between 
SAS  days  137  -  175,  a  line  from  day  175  -  207  and  another  exponential  from 
day  207  -  232.  After  this  time  the  B  channel  is  in  the  background  and  the 
hardness  ratio  is  meaningless.  Figure  4.2.3  illustrates  the  hardness  fitting 
and  figure  4.2.4  shows  confidence  contours  for  all  of  the  exponential  fits  in 
parameter  space.  The  confidence  contours  drawn  are  99%,  95%,  90%, 
68.3%  and  10%.  Table  2  tabulates  the  above  fits. 


54 


'ter  4:  Data  Analysis 


Chap 


TABLE  2a 

Decay  Fit  Parameters 

Det./Day1 

T[d-1]2 

Io[cnt/sec]3  ai[cnt/sec] 

Red.  y2 

dof 

A/137-175 

27.49 

0.11 

7768 

25 

11.4 

206 

A/200-260 

23.52 

0.05 

1666 

3.0 

8.42 

215 

B/137-175 

18.76 

0.25 

237.6 

3.0 

1.18 

226 

B/ 198-240 

16.06 

0.30 

27.16 

0.55 

1.31 

179 

H/137-175 

50.27 

1.0 

0.03736 

0.0031 

1.84 

226 

H/l 75-207 

-5.67E-5 

0.02E-5 

0.01460 

0.0 

1.06 

76 

H/207-232 

34.29 

5.3 

0.01443 

0.0011 

1.18 

140 

1.  A,  B,  and  H  refer  to  Center  Slat  A,  Center  Slat  B  and  B/A  respectively. 


Days  arc  SAS  days  =  JD  -  2442500 

2.  For  H  175-207  the  time  constant  is  the  slope  of  the  line. 

3.  For  H  175-207  Iq  is  the  y-intercept  of  the  line. 


TABLE  2b 

Ramp  &  MB  Fit  Parameters 

Day1 

m  [  d-1] 

<Jm[  d*l] 

Io[cnt/sec] 

oi[cnt/sec] 

Red.  y2 

dof 

275-285 

-2.388 

0.11 

38.86 

0.15 

4.25 

71 

285-300 

N/A 

N/A 

25.79 

0.36 

4.32 

153 

300-307 

N/A 

N/A 

32.74 

0.26 

1.37 

68 

308-338 

1.941 

0.05 

31.02 

0.36 

0.96 

248 

1. 

Days  are  SAS  days  = 

JD  -  2442500. 

Decay  Time  (days)  Decay  Time  (days)  Decay  Time  (days) 

60  62  54  18  18.6  19  19.6  -^27.2  27.4  27.6 


4.2  Results 


55 


Figure  4.2.4:  Confidence  contours  in  x2  space  for  the  various  exponential  fits.  Labels 
denote  channel  and  day  range  of  fit.  A  =  A  channel,  B  =  B  channel  and 
H  =  hardness  ratio. 
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We  now  apply  the  techniques  of  chapter  4.1  to  the  slat  collimator  and  the 
hardness  ratio.  Stellingwerf  (1978)  recommends  using  5  bins  with  2  covers 
to  uncover  periodicity.  This  binning  strategy  reduces  side  lobes  and  allows 
for  rapid  computation.  If  we  see  significant  periodicity  with  this  low 
resolution  scan  we  can  increase  the  number  of  bins  and  make  smaller  period 
steps  to  get  a  high  resolution  cut  at  the  data.  We  perform  the  analysis  for  the 
various  channels  in  figures  4.2.5a-h  using  this  strategy  in  the  ranges  7.5h  -  8h 
and  6d  -  9d. 
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Figure  4.2.5:  (a)  Stellingwerf  procedure  on  channel  A  decay  phase  at  n  =  7.75h. 

(b)  Stellingwerf  procedure  on  channel  A  decay  phase  at  FI  =  7.8d. 


Channel  A  Theta  Channel  A  Theta 
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Figure  4.2.5:  (c)  Stellingwerf  procedure  on  channel  A  ramp  phase  at  fl  =  7.75h. 

(d)  Stellingwerf  procedure  on  channel  A  ramp  phase  at  fl  =  7.8d. 
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Figure  4.2.5:  (e)  Stellingwerf  procedure  on  channel  B  decay  phase  at  n  =  7.75h. 

(f)  Stellingwerf  procedure  on  channel  B  decay  phase  at  [1  =  7.8d. 


4.2  Results 


61 


4.2.2  Monte  Carlo  Simulation 

The  values  0  obtained  from  figure  4.2.5  tend  to  support  the  hypothesis 
that  the  observed  value  of  0  is  due  to  noise.  To  obtain  an  estimate  of  the 
upper  limit  on  the  periodic  amplitude  I  conducted  a  Monte  Carlo  simulation. 

We  generated  a  thousand  random  data  sets  and  calculate  ©  at  the  period  of 
interest  for  each  set.  To  generate  the  data  sets,  for  each  data  point  we  use  the 
same  time  information  as  the  original  data  set  but  place  a  Monte  Carlo 
intensity  in  the  place  of  the  original  intensity.  We  assume  that  the  value  of  © 
obtained  from  the  real  data  is  as  high  as  it  is  because  of  noise.  The 
assumption  that  our  data  are  pure  noise  gives  us  a  worst  case  scenario  because 
we  are  trying  to  show  an  absence  of  periodicity.  Thus  any  bias  in  our 
simulation  will  favor  periodicity  and  we  err  on  the  side  of  caution.  Using 
this  assumption,  for  each  Monte  Carlo  intensity  we  generate  random 
Gaussian  noise  with  the  same  mean  and  deviation  as  the  real  data.  On  top  of 
the  noise,  we  add  a  sinusoid  with  the  amplitude  and  period  of  interest.  Figure 
4.2.6a  shows  a  sample  of  real  data  and  4.2.6b  shows  a  sample  of  simulated 
data. 

After  generating  the  random  data  sets  and  calculating  the  corresponding 
0  for  each,  we  can  make  a  histogram.  Figure  4.2.7  depicts  such  histograms 
for  the  7.75h  case  for  various  period  amplitudes,  A.  We  define  T(a)  to  be 
the  number  of  simulations  out  of  a  thousand  which  yield  a  value  of  0  as  high 
as  a  or  higher. 

To  generate  upper  limits  on  periodicity  for  a  real  data  set,  we  use  the 
value  of  ©  from  table  3  (next  page)  and  find  the  amplitude  such  that  T(0)=3 
i.e.  tne  amplitude  such  that  there  is  only  three  simulations  out  of  a 
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Figure  4.2.6:  (a)  Section  of  real  channel  A  data. 

(b)  Section  of  simulated  channel  A  data. 
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thousand  which  would  yield  that  value  or  higher.  This  corresponds  to  99.7% 
or  roughly  3a.  In  table  3,  as  explained  in  chapter  4.1,  P  is  the  probability 
that  the  observed  period  is  due  to  random  fluctuations  and  upper  limit  is  the 
percentage  periodic  amplitude  upper  limit  generated  by  the  Monte  Carlo 
simulation.  In  the  table  3a  we  choose  0  corresponding  to  the  biggest 
maximum  near  the  period  due  to  the  uncertainty  of  the  exact  value  (MB  quote 
several  periods).  We  exclude  the  improved  MB  section  of  the  data  from  table 
3a  and  will  treat  it  separately  in  chapter  4.2.3 


TABLE  3a 
Results  for  fl  =  7.8d 


Channel  Data  Section1 

e2 

P3 

3a  Upper  Limit 

A:  Decay  Phase 

.95 

.54 

2.15% 

A:  Ramp 

.955 

.67 

7.3% 

B:  Decay  Phase 

.985 

.85 

4.1% 

H:  Decay  Phase 

.987 

.87 

5.3% 

TABLE  3b 

Results  for  fl  =  7.75h 

Channel  Data  Section1 

e4 

P3 

Upper  Limit 

A:  Decay  Phase 

.995 

.95 

1.1% 

A:  Improved  MB 

.998 

.98 

6.0% 

A:  Ramp 

.998 

.98 

4.17% 

B:  Decay  Phase 

1.0 

1.0 

3.85% 

H:  Decay  Phase 

1.0 

1.0 

3.17% 

2.  From  figure  4.2.7  where  we  choose  the  biggest  max  nearby. 

3.  Two  sided  F  test  from  chapter  4.1  via  computer  (Press,  1988) 

4.  From  figure  4.2.7. 
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4.2.3  MB’s  Data 

MB  based  their  periodicity  claims  on  the  section  of  data  from  SAS  day 
285  to  322.  Figure  4.2.8a  depicts  the  Stellingwerf  procedure  carried  out  on 
this  range.  There  is  a  minimum  theta  of  .70  at  7.84d  and  the  probability  that 
this  value  of  theta  is  due  to  random  fluctuations,  P,  is  negligible.  The  percent 
amplitude,  from  figure  4.2. le  is  ~  40%.  The  region  from  day  308  to  322  is 
largely  noise.  This  section  of  data  is  basically  the  section  which  we  called  the 
ramp  and,  as  we  showed  in  section  4.2.2,  it  exhibits  minimal  periodicity.  We 
take  this  section  of  data  off  of  the  MB  data  to  get  the  improved  MB  data. 
Figures  4.2.8b  and  c  show  the  periodicity  searches  carried  out  in  these 
region.*.  We  find  a  broad  minimum  about  8.27d  with  a  theta  value  of  .575 
and  negligible  P.  This  is  the  only  section  of  data  for  which  there  is  any  strong 
periodic  behavior.  The  improved  MB  section  represents  ~22  days  or  about 
three  cycles.  We  will  discuss  the  implications  of  these  facts  in  chapter  4.3 
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Figure  4.2.8:  (a)  Stellingwerf  procedure  on  channel  A  MB  phase  at  n  =  7.8d. 

(1))  Low  res.  Stellingwerf  procedure  on  channel  A  improved  MB  phase 
(c)  High  rcs.Stellingwerf  procedure  on  channel  A  improved  MB  phase. 


4.3  Conclusions 


In  chapter  4.2  we  derived  upper  limits  on  the  periodicity  of  the  X-ray 
spectrum  of  A0620.  For  the  7.75h  period  we  found  that  the  3o  upper  limit 
on  the  amplitude  is  4.17%  in  the  A  channel  ramp.  We  showed  that  apart 
from  one  twenty  day  period,  the  3a  upper  limit  on  the  amplitude  for  the  7.8d 
period  is  7.3%. 

It  is  perhaps  not  surprising  that  we  see  no  evidence  for  the  7.75h  optical 
orbital  period.  This  period  was  determined  by  MR  through  photometric  and 
spectroscopic  observations.  There  are  no  eclipse  phenomena.  One  credible 
mechanism  which  might  impress  this  period  onto  the  X-ray  spectrum  is 
electron  scattering  from  a  high  density  electron  cloud.  If  matter  is  accreting 
via  Roche  lobe  overflow,  we  expect  such  a  cloud  near  the  LI  point  where  the 
streaming  matter  impacts  into  the  accretion  disk  (figure  4.3). 


Accretion  Disk 


Figure  4.3:  Diagram  of  electron  cloud  from  the  impact  of  the  accretion  stream  into  the 
accretion  disk  softening  photons  in  observer's  line  of  sight. 

These  electrons  would  scatter  the  X-rays  and  soften  the  resultant  spectrum. 
Checking  the  B/A  hardness  ratio  would  reveal  this  softening  to  us.  We  expect 
to  see  this  scattering  during  the  decay  phase  of  the  nova  when  mass  flow  rates 
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are  high.  If  we  saw  such  periodicity  it  would  imply  that  the  system  was 
accreting  via  Roche  lobe  overflow  and  that  we  were  at  a  fairly  high 
inclination  angle  with  respect  to  the  binary  plane  (chapter  1,  figure  1.1)  The 
fact  that  we  see  no  evidence  for  this  periodicity  suggests  that  either  the  system 
is  accreting  via  Roche  lobe  overflow,  but  is  not  at  high  inclination  angle,  or 
that  the  system  is  accreting  via  stellar  wind  capture  at  any  reasonable 
inclination  which  satisfies  the  no  eclipse  condition.  The  severe  hardness  ratio 
amplitude  3a  upper  limit  of  no  more  than  3.17%  tends  to  rule  out  high 
inclination  Roche  lobe  overflow  models. 

The  light  curves  we  obtained  are  extremely  smooth.  They  are  nearly 
perfect  exponential  decays  (Figs.  4.2.1  and  4.2.2).  Both  of  the  non-steady 
state  accretion  theories  presented  in  Chapter  1.3  predict  such  a  decay.  These 
theories  are  therefore  in  good  agreement  with  our  results. 

We  confirm  the  existence  of  three  oscillations  of  a  7.8d  period  between  7 
January,  1976  and  13  February,  1976  previously  reported  by  MB  in  the  same 
data.  This  section  of  the  data  exhibits  an  amplitude  of  approximately  40%  in 
the  A  channel.  However,  by  virtue  of  our  improved  data  reduction  based  on 
analysis  of  Crab  observations  detailed  in  chapters  2  and  3,  we  were  able  to 
show  that  in  the  rest  of  the  data  set,  in  excess  of  two  hundred  days,  the  A 
channel  has  a  3a  upper  limit  of  2.15%  before  this  section  and  after  it  the  3a 
upper  limit  in  the  A  channel  is  7.3%.  The  twenty  two  day  section  consists  of 
three  dips  and  three  rises  and  could  be  represent  3  cycles  of  a  7.8d  period.  It 
could  also  be  three  dips  and  three  rises  of  interesting,  but  periodically 
insignificant,  morphology  in  the  X-ray  light  curve.  We  therefore  find  no 
evidence  to  substantiate  the  existence  of  a  persistent  7.8d  period. 
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